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Many models of dark matter contain more than one new particle beyond those in the Standard 
Model. Often heavier particles decay into the lightest dark matter particle as the Universe evolves. 
Here we explore the possibilities that arise if one of the products in a (Heavy Particle) — > (Dark 
Matter) decay is a positron, and the lifetime is shorter than the age of the Universe. The positrons 
cool down by scattering off the cosmic microwave background and eventually annihilate when they 
fall into Galactic potential wells. The resulting 511 keV flux not only places constraints on this class 
of models but might even be consistent with that observed by the INTEGRAL satellite. 

PACS numbers: 



Although there is ample evidence for the existence of non-baryonic dark matter, the properties of the 
particle[s] that make up the dark matter are not well determined. This is not surprising as the evidence for 
dark matter to date is from observations of its gravitational effects. But this situation may change soon as 
more powerful direct and indirect detection experiments come online; indeed, there are already numerous hints 
from both sectors. Direct detection experiments DAMA/LIBRA [1], COGENT [2], CRESST [3] have events 
consistent with a dark matter signal and there are several hints of indirect detection from the Fermi Gamma 
Ray Satellite and radio observations [4, 5] . There are also two sets of observations of positrons that could be 
explained by dark matter: the observed excess of positrons over electrons in PAMELA [6] and INTEGRAL 
observations of the 511 keV line from the centre of the galaxy (see, e.g., [7-9] and [10] for a review). 

It has been suggested [11-14] that the INTEGRAL observations can be explained by dark matter decays 1 . 
The observed flux of photons ^integral — 1CU 3 ph cm~ 2 sec -1 is apparently produced by the annihilation of 
positrons with galactic electrons almost at rest. Known astrophysical sources cannot account for the totality 



of these positrons [10], so it is natural to consider positrons produced by dark matter annihilations or decays. 



In the decay scenarios considered so far, a small fraction (to/TDM) of dark matter particles decay at the 
present time into low momenta positrons (~few MeV) in the bulge, and the positrons subsequently annihilate 
with electrons to produce the 511 keV line. The required lifetime is larger than the age of the universe 
tdm ~ 10 20 sec (100 GcV /m DM ). 

Here we explore the possibility that the required positrons were produced at earlier times. For concreteness, 
we assume that the sector containing dark matter has a stable component x (which is the dark matter today) 
and an unstable component, X, that decays into x an d positrons some time after recombination. Even 
within this simple scenario, there are a number of dials to turn: the masses of x an d X; the lifetime of 
X; the charge of X; the relative abundance of X (since x is the dark matter today, its abundance is fixed 
by observations); and the branching ratios for the X decay into positrons and photons. What happens to 



1 See, however, Ref. [15] for problems with the dark matter interpretation. While the model proposed here circumvents some of 
these (the positrons being quite cold as they enter the disk and bulge), tension remains in the bulge to disk ratio. 
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the decay-produced positrons depends on the values of these parameters, and we find a rich spectrum of 
observational consequences. Thus, although motivated by the 511 keV excess, we will analyze the constraints 
and possibilities of the generic idea of a heavy particle decaying into dark matter and positrons. 

The main difference between our proposal and others is that the positrons in this class of models cool in 
the early universe very efficiently before they can annihilate. Therefore, when they enter our Galaxy, they are 
already very cold. This contrasts with models in which the positrons are produced in decays or annihilations 
of Galactic dark matter, typically concentrated near the center of the Galaxy. In those models, cooling the 
positrons before they can annihilate is a central problem. 

Section II introduces notation, the constraints, and physical processes that are relevant for any positron 
scenario. The next sections elaborate on the details of cosmological down-scattering (III); energy loss in our 
Galaxy (IV); and capture in the bulge (V). We conclude in §VI with a discussion of viable particle physics 
models. 



II. GENERAL CONSIDERATIONS AND OVERVIEW OF CONSTRAINTS 

We consider the situation where some fraction of the dark matter in the universe is unstable. The decaying 
particle is called X and it decays to \i stable dark matter, together with other decay products. We parametrize 
the number density of A-particles by its value early on (say at recombination) before decays start . Define 

nx , 



(1) 



From CMB observations, we know that the dark matter density has been approximately constant (apart from 
the expansion) since recombination. In our scenario, this is still possible if uj < 1, for then the decays produce 
little change in the dark matter density and we can safely set n x (z) = n x (z = 0)(1 + z) 3 . The number density 
of X as a function of time is therefore 

—t/r 

nx = ojn x e 1 

= 1.2 x 10- 8 cm- 3 (l + zfuj 100 GcV e-^ (2) 

m x 

where r is the lifetime of X and we have adopted the WMAP value of il c d m /i 2 = 0.11. Most of the constraints 
(and possibilities) will depend on the branching ratios, in particular what fraction of the decays produce 
positrons and photons. In this work we analyze the physics of relic positrons, therefore we will focus on 
the case where the dominant decay mode includes a positron 2 . The decay to photons is nevertheless very 
important due to the stringent constraints on diffuse gamma rays. In general, we take the branching ration to 
photons, By, as a free parameter. In some cases B 7 can be of order one. An example of this is the ncutralino 
decay into photon plus gravitino in MSSM theories with gravitino LSP. In any case, the expected value of B 7 
is at least of order a, since the Fcynman diagram of the dominant mode can be extended to have the positron 
emit a photon, leading to a decay rate a factor of a smaller. 

There is also the question of the charge of the decaying particles. The simplest case regarding the phe- 
nomcnological constraints is if X is neutral. If X is charged, the relic abundances of X + and X~ arc likely 
the same, and we call this the symmetric scenario. But there is also the possibility of an asymmetry in X + 
versus X~ (compensated by a net charge in ordinary particles so the Universe remains electrically neutral) 



2 The simplest way to ensure this is to require a small mass difference, AM < 2 GeV for neutral X or AM < 1 GeV for charged 
X. With this small mass difference, the decay to protons is kinematically forbidden and all decays produce electrons/positrons 
(with the exception of neutral pions that produce photons). A more general model with larger mass differences might also 
work but the constraints from the decay-produced protons are more complex than those we consider here. 
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similar to the observed particle anti-particle asymmetry. The constraints differ in all three of these cases 
(neutral, charged symmetric and charged antisymmetric), so we consider them separately. 
There are five sets of constraints on this scenario: 

1. Collider constraints. If X is charged and light enough, it can be produced in colliders and will leave 
a distinctive signature. 

2. Catalyzed Big Bang Nucleosynthesis (CBBN). If X is charged and the abundance is symmetric, 
the X~ would bind to light nuclei and catalyze nucleosynthesis, violating the success of the comparison 
with observed light element abundances today. 

3. Heavy Water. If X is charged both in the symmetric and asymmetric cases, the remaining X + today 
could bind with an electron and X~ with 4 He and produce heavy water, the abundance of which is 
tightly constrained. 

4. Diffuse Photon Flux from Direct Decays. A fraction B 1 of decays will produce high energy 
photons. Photons produced at z < 1000 with energies today between 100 keV and 100 GeV travel freely 
through the Universe (see, e.g., [16-18] and Fig. 4 below) so would be observed as part of the diffuse 
X-ray or Gamma-ray background today. 

5. Diffuse Photon Flux from Inverse Compton Scattering. Even if B 1 were zero, positrons (and 
electrons) produced in decays produce high energy photons via Inverse Compton Scattering off the 
Cosmic Microwave Background. These would also be part of the diffuse background today except for 
energies below 100 keV and large z where the absorption by the intergalactic medium is important. 

These constraints, and the scenario to which they apply, are summarized in Table I. 



Constraint 


Neutral X 


Charged X/Symmetric 


Charged X / Asymmetric 


Collider 


No 


Yes 


Yes 


CBBN 


No 


Yes 


No 


Heavy Water 


No 


Yes 


Yes 


Diffuse Flux: Direct Decays 


Yes 


Yes 


Yes 


Diffuse Flux: Inverse Compton 


Yes 


Yes 


Yes 



TABLE I: Constraints on three scenarios of decay-produced positrons. Asymmetric means the relic heavy X particles 
are positively charged. 

Items 1,2, and 3 have been studied carefully in the literature so we can lift results from previous work. 

Because of the presence of electromagnetic couplings, one can hope to produce these heavy long-lived 
charged particles through the well-known clcctrowcak processes. Popular examples of such particles are long- 
lived staus and gluinos i?-hadrons in the MSSM. The LHC with its two dedicated detectors is probing these 
candidates. The ATLAS collaboration [19] obtained a model- independent bound using 37 pb _1 of data at 
center-of-mass energy y/s = 7 TeV. This bound can be applied to the mass of the charged NLSP giving 
m x ± > 110 GeV. The CMS collaboration [20], with its 5 pb _1 of collected data, is able to put an even 
stronger bound on the mass mx± > 223 GeV. 

For the neutral case, for instance in e + e~ colliders like LEP, one way to look for DM is look for mono-photon 
events: e + e~ — > Efr+^. These events are produced through the same interaction responsible for decays, with 
one photon attached to e ± . More precisely the reaction is e + e~ — s- %V + 7, where x an d X appear as 
missing energy This analysis has been applied in the past to a variety of scenarios (see e.g. [21] for LEP). 
Unfortunately in our case, this strategy does not lead to a significant constraint because the operators that 
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lead to decays predict very small production rates. This conclusion holds also for hadronic machines 3 , like the 
Tevatron and the LHC where one looks for missing transverse energy accompanied by monojets pp{pp) — > E/r 
+mono-jet. See e.g. [23] for analysis on Tevatron data. 

In the charged symmetric case, the abundance at BBN should satisfy the (CBBN) constraints, which reads 

w < 2.44 x 10 3 f mx ) (3) 
VIOOGcV J \ 1CT 14 J v ' 

where Ycbbn is the maximum value allowed from catalyzed BBN for ^"-lifetimes larger than 10 5 sec. There 
are two possible values for Yqbbn available in the literature: the conservative one [24, 25] Ycbbn = 1CP 14 — 
10 -15 and the more stringent one [26, 27] Ycbbn = 10~ 16 . In both charged scenarios, charged relics 
can form heavy Hydrogen atoms that can condense in the form of anomalously heavy water in the bottom of 
oceans. In this case, the so-called heavy water bound applies to the present abundance of heavy X ± -particlcs 
[28] 

n x ± (t ) < 2.63 x 1CT 35 cm" 3 (4) 

This bound implies that the lifetime of X ± must be short enough so that, at present, almost all the X ± have 
already decayed. These bounds on charged X's are shown in the u>,t plane in Fig. 1. Note that the Heavy 
Water constraint requires lifetimes less than ~8x 10 15 s, corresponds to redshift z ~ 15. 




FIG. 1: Constraints on charged X's from Heavy Water (HW), CBBN (with the conservative upper limit Ycbbn = 
1(T 14 ) and the diffuse gamma ray background. The CBBN constraint applies only to the symmetric scenario. The 
region labelled e + e~ — > 77 would produce positrons early enough that the high background electron density would 
lead to complete annihilation and therefore no residual signal. Shaded regions are excluded, and the constraints use 
m x — 100 GeV. The horizontal blue band denotes the region in parameters space that would produce a 511 keV flux 
consistent with that observed by INTEGRAL. 



3 Notice that even though there is no coupling of dark matter to quarks, it will be generated at 1-loop level (See e.g [22]). 
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The diffuse photon flux from the decays is tightly constrained by observations, a constraint that applies 
to all scenarios. Photons can be produced directly in decays or indirectly when electrons and positrons 
scatter down off of the cosmic microwave background. We work out the details in Appendices, but here 
provide estimates for the flux due to inverse Compton scattering. (As shown in the Appendix, the constraints 
from direct decays are similar.) If the decay-produced positron has initial energy of order -Emit (typically 
of order AM/2) immediately after decay at time t = r, then the up-scattered photons will have energies 
of order 4T(r)(£ ; i n i t /m e ) 2 where T is the temperature of the CMB. These photons today will have energy 
~ 4T(r)(£'i n it/m e ) 2 /(l + z d ) = 4T Q (E init /m e ) 2 where Zd is the redshift when t = r. Therefore the diffuse flux 
will peak at 

iW-a.ekev^) 2 (5) 

independent of the lifetime of X and cut off sharply at higher energies. The diffuse flux at this peak will 
roughly be of order cux(t)/ (AttE), if these photons travel freely and are observed today. The more careful 
calculation in the Appendix, including also the attenuation through Compton scattering on cold electrons and 
photoionization, leads to 

^~^*V-4^)(^r(^f (—sec,. ,., 

Applying the SPI constraint 4 at E = 1 kcV, F < 10 cm~ 2 s _1 sr _1 keV leads to 

100GcV\ i r ^/3 

< 10 4 (t < Kr'sec) (7) 



/100GcV\ / r y 
V / V10 15 sec/ 



The full constraint (extended to larger lifetimes) is depicted in Fig. 1 when m x = 100 GeV and AM = 1 
GeV, and this holds as long as the cut-off E max is above 1 keV, or the initial positron energy is above 0.5 
GeV. The constraint loosens if Emit is smaller, so that the flux is cut-off well below a keV. 

What about the positrons themselves? The energetic positrons almost immediately inverse Compton scatter 
down off the photons in the cosmic background to become non-relativistic. This down-scattering continues 
thereafter, so early decays (zd < 20) produced positrons that thermalized with electrons and protons. These 
positrons behave as CHAMPs [29, 30] (albeit not the dominant component of the dark matter, as in the 
traditional CHAMP scenario), collapsing with baryons and dark matter to form galaxies. As we will sec 
below, they quickly annihilate in the disk or bulge of the Galaxy on timescales much too short to explain 
511 keV radiation in our Galaxy. The contribution to the 511 keV flux then is from the flux of positrons on 
highly elliptical orbits entering the disk or bulge for the first time now. In the ensuing sections, we estimate 
the 511 keV flux from these positrons, subject to the constraints depicted in Fig. 1. Our conclusion is that 
- in the regime of parameter space indicated by the blue line in Fig. 1 - positrons could produce flux with 
the same amplitude as that seen by INTEGRAL. The corollary to this is that the region above the blue band 
(uj > 10~ 5 ) is ruled out by the 511 keV line measurements. 

III. POSITRON ENERGY LOSS AFTER DECAY 

Upon production in decays, positrons/electrons lose energy mainly through scattering on CMB photons 
(Inverse Compton scattering). Taking into account expansion, energy loss is governed by [31] 

dE e+ p 2 e+ H 4_ _ ( Pe+ ^ 
dt E e 



4 http : //heasarc . gsfc .nasa. gov/docs/objects/background/dif fuse_spectrum.html 
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where ot is the Thompson cross section, H the Hubble expansion rate and p 7 is the energy density of CMB 
photons. This expression is valid as long as the center of mass energy squared, which is of order E e T, is 
smaller than m 2 or equivalcntly E e <C 10 6 GcV/(l + z), a criterion that is satisfied over the full range of 
parameter space we are considering. Using the auxiliary variable e = p e + /m e (l + z), we can write Eq. (8) in 
the more compact form 

- T (1 + zf' 2 ^e^l + zf + l with r = i ot ^° ~ 0.0113. (9) 



dln(l + z) " ' v ^ ' 3 H nU 2 



Although the full result must be obtained through numerical integration of the above equation from the 
decay redshift Zd to today, we can understand the behavior of Eq. (9) using some approximations. First, note 
that the final momentum depends only on the redshift of decay Zd, but not on the initial positron energy 
-Emit- To see this qualitatively, note that in the parameter range of interest, e starts out large, since p/m e is 
initially Ei n n/m e ~ 10 3 , and we are considering decays after z ~ 100. The scattering term in Eq. (9) is then 
of order Te(l + z) 7 ' 2 , typically very large, corresponding to rapid energy loss. In a very small redshift interval 
after production, then, e drops quickly until e(l + z) falls below one, and the square root in Eq. (9) can then 
be approximated as unity. At that point, the positron momentum is approximately equal to m e (l + z)e = m e . 
That is, inverse Compton scattering almost instantaneously slows down the positrons so that they are non- 
relativistic. 

After this steep drop, e(l + z) falls below unity, and the square root in Eq. (9) reduces to one. Integrating 
under this approximation leads to a slower but still steady decline in momentum, e cx exp{— (2/5)F[(l + 
Zdf^ 2 — (1 + z) 5 / 2 ]}. The earlier the positron is produced the more effective the loss process is, so the final 
positron momentum decreases with increasing Zd- Today, this dependence scales as e - 2T i 1 + z d) 1 /5_ 

Numerical integration exhibits these qualitative features, as shown in Fig. 2. There is the steep drop 
immediately after production, leading to positrons with momentum of order an MeV. Then the energy loss 
continues, all of which is independent of the initial energy. The earlier the positrons are produced, the 
longer the losses continue, so the final momentum is smallest for positrons produced earliest. As expected, 
positrons produced at z = 20 have momentum reduced by a factor of about 1000 more (due to the exponential 
dependence on Zd) than those produced at z = 10. Our loss equation is valid only when the positron energy is 
much greater than the CMB temperature. So, for large Zd, the final momentum we obtain from this equation 
tends to zero means that the positrons equilibrate with the photons and (p e ) ~ \JvnT . Eventually Coulomb 
scattering becomes even more important than Compton scattering and lead to the same result: low energy 
positrons that have equilibrated with the rest of the cosmic plasma. 



IV. POSITRON INTERACTIONS IN THE GALAXY 



Positrons will interact and change energy as they diffuse from the galaxy outskirts to the disk and bulge. 
As we have seen in the previous section, the relic positrons from dark matter decay thermalize immediately 
with the CMB and then after structure formation they behave as CHAMPs and are dragged to the galactic 
halo. Typical velocities for these positrons when they arrive to the disk or the bulge will be of order of 
200 km/s. Therefore the kinetic energy of these positrons will be O(0. 1) eV. Positrons and electrons with 
these energies cannot ionize or excite hydrogen atoms, so their only relevant interactions in the interstellar 
medium will be elastic Coulomb scattering with protons or free electrons and annihilation with free or atomic 
electrons. The interstellar medium (ISM) is an extremely complex environment [32] because it is made of 
regions, primarily of hydrogen, with fairly different physical properties. In Table II we present a summary 
of the different regions in the interstellar medium. While clouds of molecular (H2), atomic (HI) and ionized 
Hydrogen (HII regions) occupy a small volume fraction, the intercloud media (Hot Ionized, Warm Ionized 
and Warm Neutral) pervades most of the volume of the disk and bulge. 
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The complexity of the ISM fades away when the issue is the fate of a positron entering the disk. As we show 
here, a positron entering any of these regions will very quickly thermalize and then annihilate. To compute 
the rate of 511 keV flux from these positrons then, we can safely assume that the photons are produced very 
close to the place where the positron enters the disk/bulge. We develop this argument in this section and 
then compute the flux of positrons into the disk/bulge - and therefore the 511 keV flux - in the next section. 

Coulomb scattering between positrons and free electrons is very effective so positrons thermalize very quickly 
even in regions where the ionized fraction is relatively small. In Appendix C, we present the differential cross 
sections as function of the energy transfer and calculate the thermalization length iJtherm- The rate of positron 
energy change due to Coulomb scattering with free electrons in a medium with free electron density n e is 



Component Temperature Midplane density Filling fraction (%) Ionization fraction Scale Height 

(K) n Q (cm- 3 ) / (%) H (pc) 

Clouds 

H2 Molecular (MM) 10 - 20 

HI Cold Neutral (CNM) 50 - 100 

Traditional HII regions 8000 
Interclouds 

Warm Neutral HI (WNM) 10 3 - 10 4 

Warm Ionized HII (WIM) 8000 

Hot Ionized (HIM) ~ 10 6 



"Gaussian scale height n(z) = /no exp — (z/H) 2 . 
'Exponential scale height: n(z) = f no exp(— z/H). 

TABLE II: Typical parameters of the interstellar medium phases. 



10 2 - 10 6 0.1 < 10~ 4 75 a 

20 - 50 2 4 x 10 -4 - 10 -3 94 a 

1 - 10 5 ~ 1 30-100 6 



0.2-2 ~ 30 0.007 - 0.05 220 Q + 400 b 

0.1-0.3 ~ 20 0.6 - 0.9 900 b 

0.002 ~ 50 1 3000 a 
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Annihilation Channel Name Annihilation rate 

r ann (i) (sec -1 ) 

Charge exchange (ce) 10 -12 nn 

Radiative combination (rc) 1CP 12 n e 

Direct annihilation with free e~ (daf) 2 x 10~ 13 n e 
Direct annihilation with bound e~ (dab) 7 x 10~ 14 ng 

TABLE III: Annihilation channels for a thermal distribution of positrons with T ~ 1 eV. Values are taken from [35]. 



e + H ->• Ps + H 
e + e~ —¥ Ps 7 
e+e" -¥ 2-y 
e+H -> 27 + iT 



riven by, 

dE 



dt /3 e m e 



log = , (10) 



where (3 e is the typical electron velocity (which depends on temperature) and 6 m j n and b max are the minimum 
and maximum impact parameter. These can be related to the maximal and minimal energy transfer if max 
and -Kmin- In ionized media we can take the maximum impact parameter as 6 max = (47rn e )~ 1 / 3 ~ 1cm for 
n e ~ 1 cm~ 3 and the minimum impact parameter as 6 m i n = 2a/(m/3g) ~ 5 x 10~ 7 cm, which give a Coulomb 
logarithm of log (6 m i n /6 m ax) — 15. Dividing the typical energy (3/cbT/2) by this rate, and then multiplying 
by the velocity provides an estimate for the thermalization length 



~ 10 " 6 pc (^r) {wk) 



2 

(11) 



a scale that is much smaller than any of the regions listed in Table II. Therefore, no matter where a positron 
enters the disk or bulge, it will almost immediately thermalize with the ambient medium and then participate 
in the Galactic rotation of that cloud or inter-cloud medium. 

After thermalization, positrons will annihilate through one of the four available channels summarized in 
Table III. For each annihilation channel i with threshold energy Eth, the annihilation rate is given by 

POO 

r a „nW= / f E (E)n^ H a l {E)v(E)c\E (12) 

JE th 



where /b(-E') = 2y/ E/ir(kBT) 3 exp [— E/ksT] is the usual Maxwell-Boltzmann distribution function for 
energy and <Ji{E) and v(E) are the energy-dependent annihilation cross section and velocity respectively. 
Notice that even though the cross section for charge exchange at temperatures of T ~ cV is suppressed 
with respect to the two other channels 5 , the corresponding annihilation rate is by far the strongest. So, 
even though fewer positrons reach the threshold energy in the thermal distribution, the annihilation rate is 
still comparable to radiative combination with free electrons i.e. e + e~ — > Ps 7. The remaining channels 
occur without Ps formation and are: the well-known direct annihilation channel and annihilation with bound 
electrons that are one order of magnitude slower. Given that Ps formation channels are dominant, we can 
simply estimate the time scale over which positrons annihilate through Ps formation to be r~ nn ~ 10 12 
sec. The thermalized positrons are locked in clouds rotating around the center of the Galaxy with typical 
speed of order 200 km/sec, so the positron will annihilate before the cloud containing it has traveled 5-10 
parsec: effectively the positrons annihilate as soon as they enter the disk/bulge. From the annihilation rates 
in Table III, we can expect roughly 5-10% of the positrons to annihilate without forming Ps, depending on 



5 The threshold energy for this channel is 6.8 eV. 
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the ionization fraction of the media. This value should be compared with the latest experimental value [33] 
fp s = 0.967 ± 0.022. Furthermore, the resulting 511 keV line is broadened due to electrons' thermal velocity 
by the amount ~1.1 keV y/T/10 4 K [34], consistent with observation. 

To go further, we need to understand the extent of the disk and bulge. The interstellar medium (ISM), 
and in particular the WIM, is rather well-known close to the galactic plane [36, 37], however, it is much more 
uncertain at more than 1 or 2 kpc above the disk [38]. Therefore, in this work we will adopt a simplified 
model where we take the interstellar medium to be a cylinder of radius allowing R\ to vary between 
5 kpc < Ri < 15 kpc, and height 2h around the galactic center. Taking into account the small thermalization 
radius compared with typical size of the clouds, we do not need a precise description of the different regions 
inside this cylinder. The final ingredient in our model is the bulge which we take as an sphere of rj, ~ 1.5 kpc 
in the center of the galaxy. As we have seen, all positrons arriving to the ISM cylinder or the bulge will get 
trapped and annihilate close of the boundary of the considered region. Most of the positrons annihilating in 
the ISM cylinder will produce 511 keV photons in the disk, while in the inner 1.5 kpc, for h < 1.5 kpc, the 
produced photons will be seen as originating from the bulge. 



V. FLUX OF POSITRONS IN THE GALAXY 

As we have seen in the previous section, in practice, all the positrons entering the ISM cylinder or the 
bulge get trapped and annihilate the first time they cross the galaxy. Here we compute this incoming positron 
flux and turn it into a 511 keV flux from both the bulge and the disk. In both cases, the number density of 
positrons is governed by the rate equation: 

dn + 



dt 



-n + r ann + S , (13) 



where S is the flux of positrons into the region (either bulge or disk) and r ann is the annihilation rate. Since 
the annihilation rate is very large, the two terms on the right cancel and the number of positrons reaches 
equilibrium at 

»+ = J~- (14) 

1 ann 

To obtain the luminosity of 511 keV photons from this positron density we need to multiply by the annihilation 
rate, thereby canceling the denominator, and a factor of 2(1 — |/p s ) to account for the two photons produced 
and positronium formation. The total flux from a given region is then the volume of that region d 3 x weighted 
by 2(1 — | fp s ) S/(4:ir d 2 ) where d is the distance of the region from us. We expect then a flux of 511 keV 
photons from incoming positrons equal to 



F= 1 -* fpa [d 3 x§§,. (15) 



2tt J d(xf 

In our simplified model, the annihilations take place on the surface of the region (bulge or disk) so S has a 
delta function restricting the 3D integral over d 3 x to the 2D surface S. Other than the delta function, the 
source term is the flux, which weights the density of positrons by the the normal component of the velocity, 
so the flux of 511 keV photons reduces to 



F= tt^ 11 / dA n^ / d 3 vf+(v,x)(-v-n)e(-v-h). (16) 



Here, /+ is the occupation number of positrons, which depends on both position and velocity; n is normal to 
the surface of either the bulge or the disk; and dAa is the differential area traced out by h. The final step 
function ensures that we count only those positrons entering the disk or bulge, not those leaving. 
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The distribution of positrons in the halo is the key ingredient needed to evaluate the flux in Eq. (16). The 
simplest first estimate is to assume the positrons trace the dark matter in the halo, with the difference that 
their density is smaller by a factor of w. as defined in Eq. (2). Taking the simplest possible dark matter 
profile, an isothermal distribution, leads to 

2 

n +( r ) =UJ T-^2 — ■ ( 17 ) 
2irGr' t m x 

We fix a by requiring the local dark matter density to be equal to 0.3 GeV/cm 3 ; this translates into a = 117 
km/sec. The corresponding velocity distribution is then 

In the case of annihilating positrons, there is an additional issue to consider: positrons that have crossed 
the disk previously no longer exist as they have annihilated, so we need to count only those positrons that 
are nearing the disk or bulge for the first time now. Again we choose to model the orbits in a simple way 
to estimate this survival probability. We assume that the dark matter dominates the gravitational potential, 
which is spherically symmetric, and equal to $(r) = a 2 \n(r/r v ) with the zero-point offset r v chosen to be 200 
kpc. In this potential orbits are specified by energy E = v 2 /2 + $(r) and angular momentum L = \r x v\. 
Each orbit has two turning points (ri, r%) given by the two solutions to: 2(E — $) — L 2 /r 2 = 0, and period 

P-lT , (19) 
J ri ^2(E - $(r)) - L 2 /r 2 

The survival probability is set to one if, for given E and L, n is smaller than r (where r is the distance 
to the surface of the disk or bulge) and the period is very large; we choose the requirement P > 8 Billion 
years. Otherwise, the survival probability, s(v,x), is set to zero. Together these two constraints amount to 
the statement that only positrons on highly radial orbits are still available to fuel the 511 keV radiation today. 
It is straightforward to show numerically that these two constraints translate into separate restrictions on L 
and E. Requiring r\ < r is satisfied as long as 

L < 3.3rcr, (20) 
while the long period requirement is satisfied as long as 

E>1.2a 2 . (21) 
With these additions, the distribution function in Eq. (16) becomes 



f + (r, v = ^/2{E - $(r)), L) = f$°(r, v)s(E, L) (22) 

with 

s(E, L) = Q(E - 1.2a 2 ) e(3.3m - L). (23) 



A. Bulge 

We first compute the 511 keV flux from the bulge. For simplicity we can take all points on the surface of 
the bulge, at roughly r& = 1.5 kpc, to lie a distance of Rq = 8 kpc from us, so that d(x) comes out of the 
integral in Eq. (16). The integral over the surface has area r 2 dflft leaving 

F b = ^^tf^ jd^n J d 3 vf^(v,r b ) S (E 7 L)(-v-h)e(-v-h). (24) 
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In this spherically symmetric case, the integrand does not depend on the position on the 2D sphere, so the 
d 2 £lh integral can be done immediately giving a factor of Air. The remaining integral over the velocity is 
azimuthally symmetric so reduces to 

(1 — - fp )4irr 2 r°° 

F b = ~ -/ dvv 2 f^°{v,r b ) d6 sm6s(E,L)(-vcos6) (25) 

R® JO Jtt/2 

where the step function has been incorporated into the limits on the polar angle 9: only inward orbits, those 
with angle 9 between the velocity and the radial vector anti-aligned (cos# < 0), are included. 

To implement the constraints imposed by the survival probability, it is simplest to change dummy variables 
from (v, 9) to (E, L), where L = vr b sin9. The Jacobian is v 2 r b \ cos0|, so dividing by this, the integrals become 

F b = K £ Ps) J o dEf™(v = ^2[E-$(r b )],r b ) ^ dLL 



0.275 (3.3) 2 [1 - e- 1 ' 2 ] ^ ■ (26) 



Plugging in numbers leads to 



100GeV\ / r b 

J \1.5kpc 



F b = 217w ( — ) ( - ^— ) cm- z s- L . (27) 



This means that in order to match the observation of 10 3 cm 2 s , we would need a value of to, 

— " x Ga&v) ■ < 28 » 

B. Disk 

The flux from the disk is more complicated to compute for two reasons: (i) the distance to us now varies 
significantly as one moves along the surface of the disk, so d cannot be removed from the integral in Eq. (16) 
and (ii) the normal n is no longer parallel to the radial vector that enters the definition of the angular 
momentum. Whereas for the bulge v x r = vrsinO and v-h = vcosO, for the disk these two products involve 
different angles. Specifically, if we choose the polar axis of the d 3 v integral to lie parallel to the radial vector 
f, then L will still be equal to vrsinO, but now 

v 

v ■ h = — [R sin 6 cos <fi + h cos 6} . (29) 
r 

Here the 3D distance r = \/R 2 + h 2 where R is the cylindrical radius from the center of the Galaxy; h is the 
height of the disk; and (f> is the azimuthal angle in the velocity integration. Thus the flux coming from the 
disk is 

Fd >HW f «RR (* _ _ ,30) 



2tt Jo VR 2 + h 2 Jo R 2 + h 2 + R 2 ) - 2RR Q cosa 

x / dvv 3 f™(r,v) d9sm8s(E,L) d<j) [-£ sin cos - hcos9] Q \-Rsm9cos<t> - hcosO] . 
Jo Jo Jo 

We can simplify this equation in the limit h <C R, which is approximately true in our simplified model wherein 

the radius of the disk is much larger than its height, although we keep h in the denominator to avoid the 

singularity in R = R Q . So, 



x 



2tt Jo Jo R 2 + h 2 +R%- 2RR G cos a 

/ dvv 3 f°°(r,v) d9 sin 9 s(E,L) / d(f> [- R sin 9 cos 0} Q[-R sin 9 cos 0} . (31) 
Jo Jo Jo 
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Now 9 [— R sin 9 cost/)] selects tt/2 < (j> < 3ir/2 and the integral on <j> is trivial, becoming just 2i?sin6>. To 
proceed in this case it is more convenient to change the v integral to K = v 2 /(2a 2 ). Then, the restriction on 
L = vrsm8 < 3.3rcr can be set on K and we have, 1.2 — <f>(r)/(7 2 < K < (3.3/ sin6*) 2 /2. The upper bound 
on the K integral must be larger than the lower bound, so sin 2 6 < 3.3 2 /[2(1.2 — $(r)/<r 2 )]. For large r, $ 
is sufficiently large that this constraint is always satisfied. But for r < 3 kpc we have a bound on sin 2 9. 
Expressing the integral as one over cos 9 then leads to 

F d * 2(1 ZIM i Rd dRR ,„ J" 4*^ 



2tt J J R 2 + h 2 + R%- 2RR Q cos a (2tt) 3 / 2 

3 2 /(2sin 2 8) 

(32) 



1 /-3.3 2 /(2sin 2 8) r-a r /-3.3 2 /(2 sin 2 8) 

dcos9sin9 / K dK e~ K + I dcos9sm9 / K dK eT K 



'o r Ji.2-$(_R)/ct 2 J-l J1.2-*(fl)/ir 2 

where a r is defined as 

I r > 3 kpc 

The dimensionlcss factor inside the square brackets in Eq. (32) can be evaluated numerically; a good fit is 
C 2 (R) ^ 0.007 + 0.01( r j| H - 1) for from 1 to 20 kpc (again we neglect the small distinction here between 
r and R). Then we have: 

Fj „ ?<iziM f 4 5±<^ Cm (34) 



2 ^ ./(i? 2 + /i 2 + J R|) 2 -4 J R 2 fl| ( 27r ) 3/2 



From here, we obtain an estimate of (B/D) 1 by dividing by the flux obtained from the bulge: 

_1 = f dR - m£3 rtm (35) 

D > 2n Ji C b r 2 J{R? + R% + h?)* - 4R 2 R% n +W 



with Cb = (3.3) 2 [l — e -1 ' 2 ] ^ ~ 0.057. Integrating this expression numerically, we obtain the value of B/D. 
Given the crudeness of our disk model, we compute this for a variety of values of Rd, the radius of the disk. The 
result is shown in Figure 3. For the simple isothermal profile, the bulge to disk ratio is of order unity (in 
agreement with observations [39]) if the disk scale length is less than about 8 kpc, but drops to 0.25 for larger 
values of the disk radius. 



VI. DISCUSSION 

High energy positrons produced at early times via the decay of a second dark matter species cool down and 
will eventually get trapped in galaxies. These cooled positrons would contribute to the 511 keV flux in our 
Galaxy and - with a suitable choice of parameters - might explain the observed bulge to disk ratio of this 
flux. We conclude with a brief discussion of the possible models that can generate this kind of mass hierarchy 
and couplings in the dark sector. 

First consider the case where the unstable DM component is neutral. In supersymmetric models, while the 
possibility of two neutral species is quite natural (with both a neutralino and gravitino), the heavier species 
will often decay predominantly into photons, leading to very tight constraints from the diffuse flux [40]. This 
traces back to the fact that the lightest neutralino usually has a large photino component. More generally, we 
can use effective field theory to write down operators that lead to decays. If the dark matter species consists of 
fcrmions, the coefficient is of order A -2 where A is the UV scale above which the effective theory breaks down. 
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B/D 

2.0 r 




FIG. 3: Ratio B/D of the flux of 511 photons from the bulge and the disk for different values of the ISM disk radius 
Rd as explained in the text. 

The long lifetimes required then point to A ~ 10 11 GeV. For scalar dark matter, the operator is suppressed by 
only a single power of A, so the long lifetimes require a UV scale of order the Planck mass. In this effective field 
theory context, too, operators leading to decay to photons must be suppressed. Models that accomplish this 
typically rely on a secluded dark sector, which communicates with the Standard Model through suppressed 
interactions. In all models, it is a challenge to obtain the correct relic density since couplings are so small. 
Likely, x cannot be a conventional thermal relic and other mechanisms (such as 'freeze-in" [41] or dilution via 
a short phase of thermal inflation [42]) need to be explored. 

Getting the correct relic density is even more difficult in the case of charged dark matter. In the charge 
symmetric case, dark matter is likely to annihilate quickly into photons, leaving no appreciable X's to decay 
at late times. So a mechanism to suppress annihilations will be essential to any successful symmetric model. 
Asymmetric models do not suffer from this problem, but generating an asymmetry - which of course requires 
out-of-equilibrium - might prove challenging in a sector that interacts electromagnetically. 

Although we have touched on particle physics realizations of our scenario, in principle our constraints apply 
to any source of positrons that are operating from recombination up to now. Most of the constraints (apart 
from heavy water, colliders, and possibly direct decays to photons) apply to any mechanism that produces 
positrons in the early universe. 
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Appendix A: Diffuse Spectrum from Inverse Compton Scattering 



When high energy positrons are produced, they scatter down off of the CMB, up-scattering the photons 
to energies where they redshift freely and can be observed today as part of the diffuse background. Here we 
calculate the diffuse spectrum from this process of inverse Compton scattering. At any given time, an energy 
density Ap — E- m itnxdt/T is deposited into diffuse photons, where E- m it is the energy carried by electrons 
and positrons in a single X decay (roughly equal to tox/2 if the decay is 2-body, the branching ratio to 
positrons is unity, and mx 3> rn x ). The spectrum of photons produced by inverse Compton scattering scales 
as dp/dE oc E" 1 / 2 [43] up to a maximum energy of order E max = 4T(E in i t /m e ) 2 , so the differential spectrum 
produced by decays in a small interval dt is 



d 2 p 



-Einit^X 



dtdE 2EULe^t 



o[e b 



E] 



(Al) 



1/2 

Note that, since E max scales as -Emit, the factors of -Emit in the amplitude cancel out and the differential 
spectrum depends on Ej n i t only indirectly in the cut-off energy. This reflects the physical fact that a larger 
Emit leads to more energy injected into the diffuse spectrum per decay, but this greater energy extends to 
larger maximum energy, so that the amplitude of the spectrum at any point below the maximum does not 
depend on Ej n i t . 

To convert this to a spectrum today, use the facts that the energy today Eq is equal to E/(l + z) and the 
time interval dt is equal to — dz/[H (z)(l + z)], and integrate over all time, remembering that the left hand 
side scales as (1 + z) 3 . So, 



dp 
dEv 



1 



dz 



E- mit nx{z) 



E y 2 k H{ z ){i + zf/ 2 



9E 1/2 t 



6[E max -E (l + z)] 



(A2) 



The argument of the O function is independent of z because E max scales as 1 + z, reflecting the fact that, at 
any given time, the photons with maximum energy produced by inverse Compton scattering redshift down 
to the same maximum value of E today (4T (Ei n i t /m e ) 2 ). Dividing by the energy, we obtain the spectrum 
today: 



dn 
dEv 



m e n x (t = t) 



4T 1/2 rE 3/2 F (l + z £i ) 3 



-9 



4T (E init ) 2 



E 



Hr), 



where 



1 + z d = 1 + z(t = t) ~ 68 



10 15 sec 



2/3 



the approximate equality holding at the percent level for 1 < z < 200, and 



J(r) 



dz- 



-[t-r]/r 



(H/H )(l + z) 



9.3 x 10~ 5 f — ^ — ) 
V10 15 sec/ 
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FIG. 4: Opacity of the universe for different redshifts, (1 + z), and photon energies, Eo. The different constours 
correspond to attenuation factors, e~ r ' E °' z ^ = 0.5, 0.4, 0.3, 0.2, 0.1, 0.01 from light to dark. 

—1/2 

with the approximate equality valid for 2 < z < 600. The spectrum then scales as E with an amplitude 
governed by the dimensionless integral over redshift, which depends only on the lifetime of X. Plugging 
numbers and dividing by 47r sr leads to the present day flux in Eq. (6). 

A lingering question is whether the photons produced at large z ftravel freely through the universe or if 
they interact and are absorbed before arriving here. To answer this question we need to know the optical 
depths of the absorbtion and scattering processes. This problem was discussed in detail in Refs.[16, 44]. Using 
the code developed in [18], we find attenuation as depicted in Fig. 4. For most of the parameter space we are 
interested in here, attenuation is irrelevant, although we do include- it in the constraints shown in Fig. 1. 

Appendix B: Diffuse spectrum from radiative decays and Final State radiation 

In this appendix, we compute the diffuse spectrum of photons produced through decay of the NLSP (the 
X particles) [45]. There are two such processes that we will consider here. The first one is radiative decays 
wherein photons are produced directly. The second process is final state radiation, where a photon is emitted 
from the produced electron/positron through Brcmsstrahlung . Let us consider a decay process with partial 
decay width T 7 . Each decay leads to 7V 7 photons. For instance N 7 = 1 for X — > 7%. Following Rcf. [17], we 
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write the measured diffuse flux at redshift z as 



E^{z) = - , _ /Wi A—- ' (Bl) 



£ 7 w 47ry 2 i?(z')(l + z') (l + z') 3 /(l + z) 3 



where H{z) = Ho\/Qa + ^m(1 + z) 3 is the Hubble expansion rate and the emissivity J(Ej, z) takes the form 

J{E,z) = N^n x {z)E-^ (B2) 
where nx(-z) is the number density of the decaying particle that can be written as 

nx ( z ) = u) P^EL(l + z f e -t{*)/T 
m x 

The present photon flux from a general decay to photons can be written as 

^ - uN ±Es^ r d f e -«*)/r E ^L, (B3) 



dE 1 ' 7 4tt m x E 7 J H(z')(l + z')" dE 
where E = E(z') = (1 + z')E 1 is the energy of the photon at production. 

1. 2-body decays 

If the NLSP decays through a 2 body process X — > xii then iV 7 = 1 and dT^/dE = Br 7 TxS(E 7 (l + z') — e 7 ), 
where Br 7 = T^/Tx is the branching ratio to photons, e 7 = (l — (m x /mx) 2 ) is the photon energy at 
production and Tx = t^ 1 is the total decay rate. The delta function can be integrated using the formula 
S [f(z)] = S(z — z*)/\f'(z = where z* is the solution which satisfies f(z is ) = 0. In our case = e 7 /i? 7 — 1. 
After integration we get 



rf$ 7 _ c p c J2 cdm Br 7 e-t^O/Tjt 



F 6(e 7 - 25 7 ), (B4) 



dE 7 4?r m X Tx E^HqtJQa + O^f (e 7 /.E 7 ) 3 
where the function simply cuts energies larger than the initial energy e 7 and t(z) is given by [46] 



n 2 iog[(V» A (i + z)-z + yr> M + n A (i + z)- 3 )/^/^^] 

Taking into account that c/H ~ 1.28 x 10 28 cm , p c = 5.46 x 1(T 6 GcV/cm 3 , H Q = 72 km sec -1 Mpc -1 , 
Q\ = 0.7, Om = ri c dm + = 0.3 and f2 c( j m = 0.25, we find that the flux is 

^ = 1.12 x 10 6 cm -2 sr _1 sec -1 GcV -1 (B6) 
dE^ 

„ / GcV \ ZlOOGcV \ /l0 13 scc\ e -*(^)/^ 
xwBrJ 9(e 7 - £L) . (B7) 

This differential flux of photons must be compared to the observed diffuse gamma ray flux. The maximum of 
this flux must be smaller than that observed flux by SPI [47], COMPTEL [48] and EGRET [49]. For a given 
lifetime and AM = mx — m X i we can nno - ^ ne constraint on the product uj Br 7 as a function of the lifetime 
tx- 

As an example of 2-body decays leading to a photon, consider the decay represented in Fig. (5). Indeed, 
even though, by assumption, there is no interaction that would mediate the decay X — > \1 at tree-level, there 
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will be subdominant diagrams leading to photons 6 . The first of such diagrams is the one shown in the Figure 
(5), where one just joins the positron and electron to form a loop, where a photon is emitted. Notice that 
this decay is 2-body, so the photon will be monochromatic. The branching ratio to photons can be estimated 
as Br 7 ~ a/47r. 




FIG. 5: Diagram contributing to photon emission at 1-loop for a neutral NLSP. 



2. Final State Radiation (FSR) 



In addition to direct decays to photons, for any dark matter component decaying to a positron or a 
positron-electron pair, there are processes where photons are emitted from the final state positrons/electrons 
through Bremsstrahlung. The produced photons, which are called final state radiation (FSR) or internal 
Bremsstrahlung, will also contribute to the diffuse background of radiation. In the colinear limit, the spectrum 
of such photons at production is given by [50] 



t£T 7 
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m v , 
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(B8) 



where x = 2E/mx- In the limit of small mass splitting x = 2AM/mx "C 1, we get from Eq. (B8) 



E 



dE 



loe 



mx 



T x 
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Plugging this expression in Eq. (B3), we get the diffuse flux from internal Brehmsstrahlung 



where we have used 
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We can represent both bounds from direct decay and FSR in Fig (6); note that these bounds are consistent 
with those obtained in Refs. [51, 52]. 



Appendix C: Coulomb energy losses 

We consider here the problem of elastic Coulomb scattering of positrons on interstellar matter. The positron 
kinetic energy corresponding to a velocity of 200 km/s is of order ~ 0.1 eV so is too low to ionize or excite 



Notice that if the four-fermion interaction leading to the decay X — > \ ee contains a 75 . this process will be absent [22] . 



19 




neutral Hydrogen. Thus, at these low energies, positrons will scatter-off only free electrons or protons. Since 
the protons are very heavy, the energy transfer is completely negligible. Therefore, the energy transfer is 
dominated by Coulomb scattering of positrons and free electrons. 

Consider the elastic scattering of an electron and positron. Following Jackson [53], we work in the the 
e~ rest frame and translate results to the Lab frame. The scattering is represented in Fig. 7. The variable 
defined in the figure are b the impact parameter; V the incident angle, and 9 the scattering angle. When 
passing through its trajectory, the electron traveling at a velocity /3 will feel the electrostatic attraction given 
by Coulomb's law F = j^r. Only the force perpendicular to the trajectory will give an effect, as the parallel 
one cancels out, so the momentum transfer is 

f°° r, , a f n . , , , 2a 

Ap= F x dt = -— smipdip = —— 
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where we have changed variable using dt — fidx where da; = —bdQ/ sin 2 9. 




FIG. 7: The geometry of Coulomb scattering. 
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The energy transfer for each scattering is given by 

(A P ) 2 2c? 

AE -^T~¥^:- (C1) 

In order to obtain the total energy transfer for multiple collisions, we have to integrate over the number of 
electrons contained in the cylinder with volume element 2irbdbdx. Since dx = f3dt, we have the energy loss 
rate 

^=-AE^(2nbdb) (C2) 



where wc divided by A-k to average over all directions. Integration over the impact parameter leads to 

-nelogH^V (C3) 
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are the maximum and minimum impact parameters respectively The energy loss 
depends only logarithmically on the cutoffs so it is sufficient to obtain rough estimates of 6 max and & m i n - 
The minimum impact parameter corresponds to the closest approach and thus to the maximum energy 
transfer, so 6 m i n = -^zp- It is more difficult to define define 6 max for free electrons due the the infinite range 
of the electromagnetic force. For any given impact parameter, each scattering should take place in a volume 
AV{b) = 2irb 2 f3AT where At is the scattering duration. There should be one electron in AV(b), so the 
maximum impact parameter 6 max should satisfy n e AV r (6 max ) = I. Thus & max = (47rn e ) -1 / 3 . So, the Coulomb 
logarithm can be written as [54] 

L c = log(6 max /6 min ) ~ log(3fc B TArV2a) - 15, (C4) 

where Xd = (ksT '/ 'Aim^a) 1 / 2 is the Debye length; the screening length in a fully ionized medium, i.e. a 
plasma. 

Thus, the positron will reach the temperature of the plasma after a time tthcrm., given by, 

3k B T/2 _ 3k B T/3m e 

This is associated with a thermalization length 

i? thcrm EE pt thcrm ~ I(T 6 P C (i^) ■ ( C6 ) 



